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The origin of high-energy cosmic rays, atomic nuclei that continuously impact Earth’s atmosphere, is 
unknown. Because of deflection by interstellar magnetic fields, cosmic rays produced within the Milky 
Way arrive at Earth from random directions. However, cosmic rays interact with matter near their 
sources and during propagation, which produces high-energy neutrinos. We searched for neutrino 
emission using machine learning techniques applied to 10 years of data from the IceCube Neutrino 
Observatory. By comparing diffuse emission models to a background-only hypothesis, we identified 
neutrino emission from the Galactic plane at the 4.5c level of significance. The signal is consistent with 
diffuse emission of neutrinos from the Milky Way but could also arise from a population of unresolved 


point sources. 


he Milky Way emits radiation across the 

electromagnetic spectrum, from radio 

waves to gamma rays. Observations at 

different wavelengths provide insight into 

the structure of the Galaxy and have iden- 
tified sources of the highest-energy photons. 
For gamma rays with energies above 1 GeV, 
the plane of the Milky Way is the most prom- 
inent feature visible on the sky (Fig. 1B). Most 
of this observed gamma-ray flux consists of 
photons generated by the decay of neutral 
pions (n°), themselves produced by cosmic 
rays (high-energy charged particles) collid- 
ing with the interstellar medium within the 
Milky Way Galaxy (1). 

Photons can also be produced from inter- 
actions of energetic electrons with interstellar 
photon fields or be absorbed by ambient in- 
terstellar matter, so other messengers are 
needed to determine the cosmic-ray inter- 
actions and acceleration sites in the Galaxy. 
In addition to neutral pions, cosmic-ray inter- 
actions also produce charged pions. Charged 
pions produce neutrinos when they decay. 
Unlike photons, neutrinos rarely interact 
on their way to Earth, so they directly trace 
the location and energetics of the cosmic- 
ray interactions. 

Because both gamma rays and neutrinos 
arise from the decay of pions, a diffuse neu- 
trino flux concentrated along the Galactic 
plane has been predicted (2-5). The expected 
tera-electron volt energy neutrino flux is shown 
in Fig. 1D, calculated from the giga-electron 
volt gamma-ray observation (7). In addition to 
the predicted diffuse emission, the Milky Way 
is densely populated with numerous high- 


*Corresponding author: analysis@icecube.wisc.edu; F. Halzen 
(francis.halzen@icecube.wisc.edu) 

tlceCube Collaboration authors and affiliations are listed in the 
supplementary materials. 


IceCube Collaboration, Science 380, 1338-1343 (2023) 


energy gamma-ray point sources (also visible 
in Fig. 1B), several classes of which are po- 
tential cosmic-ray accelerators and therefore 
possible neutrino sources (6-10). This makes 
the Galactic plane an expected location of 
neutrino emission. 

Previous searches for this signal by using 
neutrino detectors (71-14) have not found a 
statistically significant signal (P values 20.02). 
The development of deep learning techniques 
in data science has produced tools (15-17) that 
can identify a larger number of neutrino inter- 
actions in detector data, with improved angu- 
lar resolution over earlier methods. We applied 
these deep learning tools to data from the 
IceCube Neutrino Observatory to search for 
evidence of neutrino emission from the Galac- 
tic plane. 


Cascade events in IceCube 


The IceCube Neutrino Observatory (18), located 
at the South Pole, is designed to detect high- 
energy (21 TeV) astrophysical neutrinos and 
identify their sources. The detector construc- 
tion, which deployed instruments within a 
cubic kilometer of clear ice, was completed 
in 2011; 5160 digital optical modules (DOMs) 
were placed at depths between 1.5 and 2.5 km 
below the surface of the Antarctic glacial ice 
sheet. Neutrinos are detected through Cheren- 
kov radiation, emitted by charged secondary 
particles that are produced by neutrino inter- 
actions with nuclei in the ice or bedrock. Be- 
cause of the large momentum transfer from 
the incoming neutrino, the directions of sec- 
ondary particles are closely aligned with the 
incoming neutrino direction, enabling the iden- 
tification of the neutrino’s origin. 

The two main detection channels are cas- 
cade and track events. Cascades are short- 
ranged particle showers, predominantly from 
interactions of electron neutrino (v,) and tau 
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muon neutrino (v,,), and v,] on nuclei. Be- 
cause the charged particles in cascade events 
travel only a few meters, these energy deposi- 
tions appear almost point-like to IceCube’s 
125-m (horizontal) and 7- to 17-m (vertical) 
instrument spacing. This results in larger di- 
rectional uncertainties than tracks. Tracks are 
elongated energy depositions (often several 
kilometers long), which arise predominantly 
from muons generated in cosmic-ray particle 
interactions in the atmosphere or muons pro- 
duced by interactions of v, with nuclei. The 
energy deposited by cascades is often con- 
tained within the instrumented volume (un- 
like tracks), which provides a more complete 
measure of the neutrino energy (19). 

Searches for astrophysical neutrino sources 
are affected by an overwhelming background 
of muons and neutrinos produced by cosmic- 
ray interactions with Earth’s atmosphere. At- 
mospheric muons dominate this background; 
IceCube records about 100 million muons for 
every observed astrophysical neutrino. Whereas 
muons from the Southern Hemisphere (above 
IceCube) can penetrate several kilometers deep 
into the ice, muons from the Northern Hemi- 
sphere (below IceCube) are absorbed during 
passage through Earth. Because of this shield- 
ing effect, and the superior angular resolution 
of tracks over cascades (<1° compared with 
<10°, respectively; both above 10 TeV), searches 
for neutrino sources that use IceCube typically 
rely on track selections, making them most 
sensitive to astrophysical sources in the North- 
ern sky (20). 

The Galactic Center, as well as the bulk 
of the neutrino emission expected from the 
Galactic plane, is located in the Southern sky 
(Fig. 1, C and D). To overcome the muon back- 
ground in the Southern sky, analyses of IceCube 
data often use events in which the neutrino 
interaction is observed within the detector 
(21, 22). The selection of these events greatly 
reduces the background rate of cosmic-ray 
muons that enter the instrumented volume 
from the Southern sky. Unlike these atmo- 
spheric muons, atmospheric neutrinos (23) 
generally cannot be distinguished (by their in- 
teractions in the detector) from astrophysical 
neutrinos. Nevertheless, with increasing ener- 
gy, an increasing fraction of the atmospheric 
neutrinos from the Southern sky (above IceCube) 
can be excluded by eliminating events with 
simultaneous muons that originate from the 
same cosmic-ray air-shower that produced the 
atmospheric neutrino (24, 25). At the tera- 
electron volt energies relevant for searches of 
Galactic neutrino emission, the majority of 
these atmospheric neutrinos remain as a sub- 
stantial background in searches for astrophys- 
ical neutrinos. This background is dominated 
by muon neutrinos, which are largely detected 


1 of 6 


RESEARCH | RESEARCH ARTICLE 


m 
u 


y Optical 


Latitude [b] 
oO 


A 

15° 
= 15 20> 
Q y Gamma Ray. 2 
v a 
Z 0 Š 
2 Ss 
5 a 
15° 0.52 
= 15 7 4 
z Predicted r° Northern Sky / 
oS / 2T 
EE yaen i 205 
È / \ z7 
5 / Southern Sky Southern Sky \ 

-15° 1 : 0 
= ib 0.5 — 
Q v Analysis Expectation | 
v 5 
S 0 Typical Event Uncertainty £ 
2 i © 
EEN D : 

-15° 0.0 = 

ES” 4 


ae 
— a 


$ 


Latitude [b] 
e] 


180° 120° 60° 0° 


Galactic Coord. 


N 
Significance 
[o] 


Oo 


-60° -120° -180° 


Galactic Longitude [/] 


Fig. 1. The plane of the Milky Way Galaxy in photons and neutrinos. (A) to 
(E) are in Galactic coordinates, with the origin being at the Galactic Center, 
extending +15° in latitude and +180° in longitude. (A) Optical color image (39), 
which is partly obscured by clouds of gas and dust that absorb optical photons. 
[Credit: A. Mellinger, used with permission.] (B) The integrated flux in gamma 
rays from the Fermi Large Area Telescope (Fermi-LAT) 12-year survey (40) 

at energies greater than 1 GeV, obtained from the Fermi Science Support Center 
and processed with the Fermi-LAT ScienceTools. (C) The emission template 
calculated for the expected neutrino flux, derived from the n° template that 


matches the Fermi-LAT observations of the diffuse gamma-ray emission (1). 
(D) The emission template from (C), after including the detector sensitivity to 
cascade-like neutrino events and the angular uncertainty of a typical signal event 
(7°, indicated by the dotted white circle). Contours indicate the central regions 
that contain 20 and 50% of the predicted diffuse neutrino emission signal. 

(E) The pretrial significance of the IceCube neutrino observations, calculated 
from the all-sky scan for point-like sources by using the cascade neutrino event 
sample. Contours are the same as in (D). Gray lines in (C) to (E) indicate the 
northern-southern sky horizon at the IceCube detector. 


as tracks in IceCube. The selection of cascade 
events instead of track events therefore reduces 
the contamination of atmospheric neutrinos— 
by about an order of magnitude at tera-electron 
volt energies—and permits the energy thresh- 
old of the analysis to be lowered to about 1 TeV. 

In the Southern sky, the lower background, 
better energy resolution, and lower energy 
threshold of cascade events compensate for 
their inferior angular resolution, compared 
with those of tracks. This is particularly true for 
searches for emission from extended objects, 
such as the Galactic plane, for which the size 
of the emitting region is larger than (or similar 
to) the angular resolution. Compared with 
track-based searches, cascade-based analyses 
are more reliant on the signal purity and less 
on the angular resolution of individual events. 
We therefore expect analyses based on cascades 
to have substantially better sensitivity to ex- 
tended neutrino emission in the tera-electron 
volt energy range from the Southern sky. 


Application of deep learning to cascade events 


To identify and reconstruct cascade events in 
IceCube, we used tools based on deep learn- 
ing. These tools are designed to reject the 
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overwhelming background from atmospheric 
muon events, then to identify the energies and 
directions of the neutrinos that generated the 
cascade events. IceCube observes events at a 
rate of about about 2.7 KHz (18), arising mostly 
from background events (atmospheric muons 
and atmospheric neutrinos) that outnumber 
signal events (astrophysical neutrinos) at a 
ratio of roughly 10°:1. To search for neutrino 
sources, event selection was required to im- 
prove the signal purity by orders of magnitude. 

Previously used event selections for cascade 
events (22, 26, 27) relied on high-level observ- 
ables, such as the event location within the 
IceCube volume and total measured light levels, 
to reduce the initial data rate. In subsequent 
selection steps, more computing-intensive se- 
lection strategies were performed, such as the 
definition of veto regions within the detector, 
to further reject events identified as incoming 
muons. We adopted a different approach, 
using tools based on convolutional neural net- 
works (CNNs) (15, 28) to perform event selec- 
tions. The high inference speed of the neural 
networks (milliseconds per event) allowed us 
to use a more complex filtering strategy at 
earlier stages of the event selection pipeline. 
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This retains more low-energy astrophysical 
neutrino events (Fig. 2) and includes cascade 
events that are difficult to reconstruct and dis- 
tinguish from background because of their lo- 
cation at the boundaries of the instrumented 
volume or in regions of the ice with degraded 
optical clarity (from higher concentrations of 
impurities in the ice). 

After the selection of events, we refined 
event properties, such as the direction of the 
incoming neutrino and deposited energy, using 
the patterns of deposited light in the detector. 
The likelihood of the observed light pattern 
under a given event hypothesis was maximized 
to determine the event properties that best 
describe the data. For this purpose, we used 
a hybrid reconstruction method (J6, 17) that 
combines a maximum likelihood estimation 
with deep learning. In this approach, we used 
a neural network (NN) to parameterize the 
relationship between the event hypothesis 
and expected light yield in the detector. This 
smoothly approximates a (more computation- 
ally expensive) Monte Carlo simulation while 
avoiding the simplifications that limit other 
reconstruction methods (19, 29). Starting with 
an event hypothesis, the NN models the photon 
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Fig. 2. Neutrino effective area and event selection comparison. (A) The all-flavor southern sky effective area 
(Aer) of the IceCube dataset, averaged over a solid angle in the declination (5) range between -90° and -5° 

as a function of E, the true neutrino energy. Results are shown for the deep learning event selection used in this 
work (dark blue), a previous cascade event selection (light blue) (12), and a previous track event selection (gray) 
(20) applied to the IceCube data. (B) The number of expected signal events (Nastro) in the Southern sky per 
energy bin per year for each event selection, assuming an isotropic astrophysical flux (22). Calculations are based 
on equal contributions of each neutrino flavor at Earth because of neutrino oscillations. 


yield at each DOM. Symmetries (such as rota- 
tion, translation, and time invariance of the 
neutrino interaction) and detector-specific do- 
main knowledge are exploited by directly in- 
cluding them in the network architecture, which 
is analogous to how a Monte Carlo simulation 
would exploit this information. This differs 
from previous CNN-based methods used in 
neutrino telescopes (15), which inferred the 
event properties directly from the observed 
data. However, the observed IceCube data 
are already convolved with detector effects, 
making it difficult to exploit the underlying 
symmetries. Our hybrid method is intended 
to provide a more complete use of available 
information. A description of the hybrid meth- 
od has been published previously (16), and 
we discuss its application to our dataset (30). 

We found that this deep learning event se- 
lection retains more than 20 times as many 
events as that retained with the selection meth- 
od used in the previous cascade-based Galactic 
plane analysis of IceCube data (Fig. 2) (72). It 
also provides improved angular resolution, by 
up to a factor of 2 at tera-electron volt energies 
(fig. S5) (16). The increased event rate is mostly 
due to the reduced energy threshold and the 
inclusion of events near the boundaries of 
the instrumented volume (fig. S3). We analy- 
zed 10 years of IceCube data, collected be- 
tween May 2011 and May 2021. A total of 
59,592 events were selected over the entire 
sky in the energy range of 500 GeV to several 
peta-electron volts, compared with 1980 events 
from 7 years in the previous selection (12). We 
estimate that the remaining sample has an 
atmospheric muon contamination of about 6% 
(30), whereas the astrophysical neutrino con- 
tribution is estimated to about 7%, assuming 
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the observed flux (22). The remaining 87% of 
the events are atmospheric neutrinos. These 
fractions are not used in the analysis directly; 
instead, we used the entire sample to derive a 
data-driven background estimate. 


Searches for Galactic neutrino emission 


We used this event selection to perform 
searches based on several neutrino emission 
hypotheses (30). For each hypothesis, we used 
a previously described maximum likelihood- 
based method (31), modified to account for 
signal contamination in the data-derived back- 
ground model (77, 12). These techniques, decided 
a priori and blind to the reconstructed event 
directions, infer the background from the data 
itself, avoiding the uncertainties introduced by 
background modeling. We calculated P values 
by comparing the experimental results with 
mock experiments performed on randomized 
experimental data. The backgrounds for these 
searches—consisting of atmospheric muons, 
atmospheric neutrinos, and the flux of ex- 
tragalactic astrophysical neutrinos—are each 
largely isotropic. The rotation of Earth ensures 
that for a detector located at the South Pole, 
the detector sensitivity to neutrinos at differ- 
ent right ascensions is fairly uniform in each 
declination band. Therefore, we estimated back- 
grounds by scrambling the right ascension 
value of each event, preserving all detector- 
specific artifacts in the data. Any systematic 
differences between the modeling of signal 
hypotheses and the true signal could reduce 
the sensitivity of our search but would not 
invalidate the resulting P values. 

The source hypothesis tests were defined a 
priori. They include tests for the diffuse emis- 
sion expected from cosmic rays interacting 
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with the interstellar medium in the Galactic 
plane, tests that use catalogs of known Galac- 
tic sources of tera-electron volt gamma rays, 
and a test for neutrino emission from the 
Fermi Bubbles (large areas of diffuse gamma- 
ray emission observed above and below the 
Galactic Center) (32). We also performed an all- 
sky point-like source search and a test for emis- 
sion from a catalog of known giga-electron volt 
(mostly extragalactic) gamma-ray emitters (sup- 
plementary text). The results for each test (30) 
are summarized in Table 1. 


Galactic plane neutrino searches 


We tested three models of Galactic diffuse 
neutrino emission, extrapolated from the ob- 
servations in gamma rays (Fig. 1B). These mod- 
els are referred to as n°, KRA®, and KRA®? (33) 
and are each derived from the same under- 
lying gamma-ray observations (7). The model 
predictions depend on the distribution and 
emission spectrum of cosmic-ray sources in 
the Galaxy, the properties of cosmic-ray diffu- 
sion in the interstellar medium, and the spa- 
tial distribution of target gas. Each neutrino 
emission model was converted to a spatial tem- 
plate, then convolved with the detector ac- 
ceptance and the event’s estimated angular 
uncertainty, to produce an event-specific spatial 
probability density function (shown for a typical 
event angular uncertainty of 7° in Fig. 1D). 

The x° model assumes that the mega-electron 
volt-to-giga-electron volt n? component, infer- 
red from the gamma-ray emission, follows a 
power law in photon energy (E) of E” and 
can be extrapolated to tera-electron volt en- 
ergies with the same spatial emission profile. 
The KRA, models include a variable spectrum 
in different spatial regions, use a harder (on 
average) neutrino spectrum than that of the x° 
model, and include a spectral cutoff at the 
highest energies (33). In this analysis, the KRA, 
models are tested with a template that uses a 
constant, model-averaged spectrum over the 
sky, roughly corresponding to an E °” power 
law, with either a 5 or 50 PeV cosmic-ray en- 
ergy cutoff for the KRA® and 50 models, 
respectively. The KRA, models predict more 
concentrated neutrino emission from the Ga- 
lactic Center region, whereas the n°? model 
predicts events more evenly distributed along 
the Galactic plane. The corresponding neutrino 
spectrum predicted by each of these models has 
a cutoff at about 10 times lower energies. 

We performed Galactic template searches 
with the same methods as those of previous 
Galactic diffuse emission searches (J1, 12). 
Because of the uncertainties in the expected 
distribution of sources, and their emission spec- 
trum and cosmic-ray diffusion, we make no 
assumption about the absolute model nor- 
malization. Instead, the analyses include an 
unconstrained free parameter for the number 
of signal events (n,) in the entire sky, which 
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Fig. 3. Galactic plane test-statistic contributions. The contribution to the test-statistic t is shown in galactic coordinates (longitude / and latitude b) for each of 
the three tested Galactic plane models. The overall test-statistic value was obtained through integration over the whole sky. (A to C) The contributions of each 
model for the observed data. (D to F) The contributions of each model for a single randomly selected mock experiment by using scrambled data. In (A) to (F), 
contours enclose 20% (white) and 50% (gray) of the predicted model flux; for the x° model, these are the same as in Fig. 1, D and E. The 50% contours contain about 
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provides flux normalization, while keeping the 
spectrum fixed to the model. Results for each 
model are summarized in Table 1. We rejected the 
background-only hypothesis with significances 
of 4.710, 4.376, and 3.960 for the x°, KRAŽ, and 
50 models, respectively. Although these 
three hypotheses are correlated, we applied a 
conservative trial factor of 3 to the most signif- 
icant of these values, leading to a trial-corrected 
P value equivalent to a significance of 4.480. 
The best-fitting fluxes are also listed in Table 1. 
The flux normalization of the n? model is 
quoted at 100 TeV, assuming a single power 
law; however, the KRA, models have a more 
complex spectral prediction and are therefore 
quoted as multiples of the predicted model 
flux. These fluxes correspond to best-fitting 
values of 748, 276, and 211 signal events (ns) 
in the IceCube dataset for the n°, KRAŠ, and 
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KRA? models, respectively. A visualization 
of the template results is shown in Fig. 3, A 
to C, in the form of a map of the per-steradian 
contribution to the results in the sky region 
near the Galactic Center for each of the Ga- 
lactic plane models. Similar maps for a ran- 
domly selected mock experiment are also 
shown for comparison (Fig. 3, D to F). 

An all-sky point source search was also per- 
formed (supplementary text), in which the sky 
was divided into a grid of equal solid-angle 
bins, spaced 0.45° apart, and each point was 
tested as a neutrino point source. The resulting 
significances are shown in Fig. 4. Some loca- 
tions have excess emission over the background 
expectations, including some in spatial coinci- 
dence with known gamma-ray emitters, such as 
the Crab Nebula, 3C 454.3, and the Cygnus X 
region. However, after accounting for trial fac- 


30 June 2023 


Fig. 4. All-sky point source search. A map of 

the best-fitting pretrial significance for the all-sky 
search, shown on an Aitoff projection of the celestial 
sphere, in equatorial coordinates (J2000 equinox). 
The Galactic plane is indicated with a grey curve, 
and the Galactic Center is indicated with a gray 
dot. Although some locations appear to have 
significant emission, the trial factor for the number 
of points searched means that these points are 

all individually statistically consistent with 
background fluctuations. 


tors, no single point in the map is statistically 
significant (Table 1). This also implies that the 
emission that is present in the Galactic template 
analyses is not due to a single point source. 


Searches using catalogs of Galactic sources 


The total gamma-ray signal from the Galactic 
plane includes a contribution from several 
strong gamma-ray point sources (1). We there- 
fore searched for correlated neutrino emission 
from three distinct catalogs of Galactic sources. 
Previous work had classified each source as a 
supernova remnant (SNR), pulsar wind nebula 
(PWN), or other unidentified (UNID) Galactic 
source, based on observations in tera-electron 
volt gamma rays (34, 35). Under the assump- 
tion that multiple sources in each class emit 
neutrinos, stacking these sources in a single 
analysis provides higher sensitivity compared 
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Fig. 5. Energy spectra for 
each of the Galactic plane 
models. Energy-scaled, sky- 
integrated, per-flavor neutrino 
flux is shown as a function of 
neutrino energy (E,) for each of 
the Galactic plane models. 
Dotted lines are the predicted 
values for the x° (dark blue), 
KRA® (orange), and KRA-” (light 
blue) models. Solid lines are our 
best-fitting flux normalizations 
from the IceCube data. Shaded 
egions indicate the lo uncer- 
tainties; they extend over the 
energy range that contributes 
to 90% of the significance. 
These results are based on the 
all-sky (4r sr) template and are 
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Table 1. Summarized results of the neutrino emission searches. The flux sensitivity and best-fitting 
flux normalization (®) are given in units of model flux (MF) for the KRA, templates and for the n? analyses 
as [E aN at 100 TeV, in units of 107? TeV cm? s (where an is the differential number of neutrinos per 
flavor, N, and neutrino energy, E). P values and significances are calculated with respect to the 
background-only hypothesis. Pretrial P values for each individual result are listed for the three diffuse 
Galactic plane analyses and three stacking analyses, and posttrial P values are given for the other analyses 
(supplementary text). Because of the spatial overlap of the stacking catalogs with the diffuse Galactic 
plane templates, strong correlations between these searches are expected. More detailed results for each 


search are provided in tables S1 to S5. 


Flux sensitivity ® 


P value Best-fitting flux ® 


Diffuse Galactic plane analysis 


0.11 x MF 


0.55408 x MF 


3.72 x 10° (3.960) 0.377 on x MF 


5.90 x 10-4 (3.246)* 


3.39 x 10-4 (3.406)* 


Hotspot (south) 


0.06 (1.520) 


0.28 (0.580) 


0.46 (0.100) 


*Significance values that are consistent with the diffuse Galactic plane template search results. 


with individual source searches, because the 
neutrino fluxes add together, whereas random 
background adds incoherently (36). The ob- 
jects in each catalog were selected according 
to the observed gamma-ray emission above 
100 GeV and the detector sensitivity, following 
previously described methods (20). We chose 
the 12 sources from each category with the 
strongest expected neutrino flux and weighted 
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them under the hypothesis that each contrib- 
utes equally to the flux (supplementary text). 
The total number of signal events and the 
spectral index are left as free parameters for 
each catalog search. The resulting P value for 
each catalog search is shown in Table 1. Each 
result rejects the background-only hypothesis 
at the 30 level or above. However, we do not 
interpret these neutrino event excesses as a 
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detection because the objects in these Galactic 
source catalogs overlap spatially with regions 
that predict the largest neutrino fluxes in the 
Galactic plane diffuse emission searches. 


Implications of Galactic neutrinos 


The neutrino flux we observed from the Galac- 
tic plane could arise from several different 
emission mechanisms. The predicted energy 
spectra integrated over the entire sky is shown 
in Fig. 5 for each of the Galactic plane models 
and their best-fitting flux normalization. Model- 
to-model flux comparisons depend on the 
regions of the sky considered. The KRA, best- 
fitting flux normalizations are lower than pre- 
dicted, which could indicate a spectral cutoff 
that is inconsistent with the 5 and 50 PeV 
values assumed. The simpler extrapolation of 
the x° model from giga-electron volt energies 
to 100 TeV predicts a neutrino flux that is a 
factor of ~5 below our best-fitting flux. How- 
ever, the best-fitting flux for the z? model ap- 
pear to be consistent with recent observations 
of 100-TeV gamma rays by the Tibet Air Shower 
Array (fig. S8) (37). The x° model mismatch 
could arise from propagation or spectral differ- 
ences for cosmic rays in the Galactic Center 
region, or from contributions from unresolved 
neutrino sources. 

We used model injection tests to quantify 
the ambiguity between different source hy- 
potheses. In these tests, the best-fitting neu- 
trino signal from one source search was 
simulated, then the expected results in all 
other analyses were examined. Injecting a 
signal from the x° model analysis, with a flux 
normalization equal to the best-fitting value 
from the observations, produces a median sig- 
nificance that is consistent with the best-fitting 
values for all other tested hypotheses (within 
the expected statistical fluctuations). This in- 
cludes the 30 excess observed in Galactic source 
catalog searches. Individually injecting the 
best-fitting flux of any one of the tested Ga- 
lactic source catalogs, at the flux level observed, 
did not recover the observed n° or KRA, model 
results. However, the angular resolution of the 
sample and the small number of equally 
weighted sources included in these catalogs 
does not constrain emissions from these broad 
source populations. It is plausible that many 
independently contributing sources from the 
Galactic plane could show a similar result to 
diffuse emission from interactions in the inter- 
stellar medium. These tests favor a neutrino 
signal from Galactic plane diffuse emission, 
but we do not have sufficient statistical power 
to differentiate between the tested emission 
models or identify embedded point sources. 

The neutrinos observed from the Galactic 
plane contribute to the all-sky astrophysical 
diffuse flux previously observed by IceCube 
(Fig. 5) (21, 22, 38). The fluxes we infer for each 
of the Galactic template models contribute 
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~6 to 13% of the astrophysical flux at 30 TeV. 
These comparisons are complicated because 
of different spectral assumptions and tested 
energy ranges used in each analysis. Addition- 
ally, the observed Galactic flux is integrated 
over the entire sky, so local emission along the 
central region of the Galactic Plane contrib- 
utes a higher fraction to the total flux. 

The observed excess of neutrinos from the 
Galactic plane provides strong evidence that 
the Milky Way is a source of high-energy neu- 
trinos. This evidence is consistent with the 
expected interactions of cosmic rays within 
the Galaxy, as established with gamma-ray mea- 
surements. It complements IceCube’s measure- 
ments of the diffuse extragalactic flux to provide 
a more complete picture of the neutrino sky. 
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Editor’s summary 

Observations of high-energy astrophysical neutrinos have shown that they mostly originate from extragalactic sources 

such as active galaxies. However, gamma ray observations show bright emission from within the Milky Way galaxy, 
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and astrophysical gamma rays and neutrinos are expected to be produced by the same physical processes. The 
IceCube Collaboration searched for neutrino emission from within the Milky Way (see the Perspective by Fusco) and 
found evidence of extra neutrinos emitted along the plane of the Galaxy, which is consistent with the distribution of 
gamma-ray emission. These results imply that high-energy neutrinos can be generated by nearby sources within the 
Milky Way. —Keith T. Smith 


View the article online 
https:/Awww.science.org/doi/10.1126/science.adc9818 
Permissions 
https://www.science.org/help/reprints-and-permissions 


Use of this article is subject to the Terms of service 


Science (ISSN ) is published by the American Association for the Advancement of Science. 1200 New York Avenue NW, Washington, DC 
20005. The title Science is a registered trademark of AAAS. 

Copyright © 2023 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim 
to original U.S. Government Works 


